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Abstract To achieve a high precision 7 mass measurement at BESIII, two free para- 
meters (m, and €) and three parameters (m,, € and opga) fittings are simulated using 
sampling technique. For two parameters fitting, two points should be taken, the one is 
near the threshold of t*7~ production to determine m+, the other point is at 3.595 GeV 
to determine efficiency. The ratio of luminosity at the two points is 3 to 1. For three 
parameters fitting, one point should be added at the low energy region with about 10% of 
the total luminosity. The optimal ratio of luminosity at the other two points is still 3 to 


1. 
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1 Introduction 


The mass of 7 lepton is a fundamental parameter in the standard model, measurement 
of Tr mass can be used to provide a significant test of lepton universality. Pseudomass 
technique [1, 2] and threshold scan method are employed to measure the r mass. The 
former relies on the reconstruction of the invariant mass and energy of the hadronic 
system in hadronic 7 decays, the latter needs a good understanding of the production 
cross section near threshold. 

A most precise measurement of Tm has been performed by BES collaboration [3, 4, 5] 


using threshold scan method: 


mM, = 1776961) 5 toe MeV/c’, (1) 
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with a relative error of 1074, which dominates the PDG 7 mass value. 

Recently, Some theoretical calculations have achieved the accuracy at the level of 1074 
[6, 7] for the production cross section near threshold. Experimentally, the high precision 
of beam energy, about 107° GeV [8], has been determined by KEDR collaboration using 
depolarization technique. Large T data sample is expected at the forthcoming experiment 
BESIII [9]. Therefore, we have great interest to wonder what precision of r mass will be 
obtained in the near future. 

In the following parts, we begin with an introduction of methodology, then present 
one parameter optimization conclusion. The two parameters and three parameters opti- 
mization are researched in detailed. At last, possible effect of luminosity and time on a 


certain statistic error of r mass is discussed. 


2 Methodology 


Within a specified period of data taking time or equivalently for a given integrated 
luminosity, we try to find out the scheme which can provide the highest precision on m,. 
The sampling technique is utilized to simulate various data taking possibilities among 
which the optimal one is chosen. For certain simulation, the likelihood function is con- 
structed as follows [3, 4, 5]: 


LF(m,) = [[ = —., (2) 


where N; is the observed number of t+7~ events for ey final state! at scan point i; ju; is 


the expected number of events and given by 
pilm) = [e - Bep ` Tobs(Mr, Lim) + OBa] ` Li - (3) 


In Eq. 3, £; is the integrated luminosity at the i-th point; € is the overall efficiency of 
eu final state for identifying 7*7~ events, which includes trigger efficiency and event 
selection efficiency; Beu is the combined branching ratio for decays T* — ety,V, and 
T — W D,V,, or the corresponding charge conjugate mode; Cobs, which can be calculated 
by the improved Voloshin’s formulas |6], is the observed cross section measured at energy 
Et, with m, as a parameter; and ogg is the total cross section of background channels 
after 7*7~ selection. 


As a tentative beginning, the energy interval to be studied is divided evenly, viz. 
B,= Eo+ (i—1)x ôE, (¢=1,2,..., Np) (4) 
where the initial point Eo = 3.50 GeV, the final point Ey = 3.595 GeV, and the fixed step 
OF = (Ey — Eo) /Np with Np being the number of energy points. For a given integrated 


luminosity (Liot) it is also apportioned averagely at each point, Li = Ltot/ Not. 


'For briefness, the ep channel means Tt > etveDr, T7 —> LW DuVr, OF TT TV, T > e Vez . 
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Input FE; and £; 
sampling by p;(7) 


Repeat fitting 


Calculate mŻ and S$, 


Figure 1: Flow chart of sampling simulation, where 7 (i = 1,2,---, Np) indicates certain 


scheme and j (i = 1,2,---, Nsamp) sampling times. 


In order to reduce the statistical fluctuation, sampling is repeated many times (Nsamp = 
200 for our study) for each scheme (say for each N,,), the average value and corresponding 


variance of the fit out 7 mass are worked out as follows [11] : 


, { “eer 
m ae a (5) 
A = + E (6) 
a= aT È ons — 


Here it should be noted that 7 indicates the certain scheme, whose value can be 1 while 7 
indicates the sampling times which equals to 200 in the following study. Without special 
declaration, the meaning of the average defined by Eqs. (5) and (6) will be kept in the 
study follows. 

The general flow chart of sampling and fitting research is presented in Fig. 1. 


During the study of high precision of m, mass measurement, we mainly concern with: 
1. the number of points to be taken; 
2. the optimal position of the data points; 


3. the relation between the luminosity and the precision. 


3 Conclusion of one-parameter fitting 


In one parameter fitting, m, is the solo free parameter. The region that the derivative 
of cross section falls to 75% of its maximum value, is proved to be effective to decrease 
the error of r mass. The total luminosity of 45 pb~! is rationed averagely in the sensitive 
region from one to six points. The simulation indicates that one point suffices to give rise 
to small uncertainty. The sensitive region is scanned to hunt for the optimal position, the 
smallest error of r mass is obtained near the m, threshold, at 3.55379 GeV. One point 
in the optimal location with luminosity of 63pb~! is sufficient for a statistical accuracy 
better than 0.1 MeV. 


4 ‘'Two-parameter fitting 
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Figure 2: The variation of Ae and S, with the number of points Np. The circles and stars 
denote the values of S, and Ae respectively. The empty legends denote the the luminosity 
allotted evenly scheme, and the filled legends denote the expected event number equal 


scheme. 


For two parameters fitting, m, and e are free variables. The optimal point, 3.55379 
GeV, which determines the r mass most economically in one parameter fitting [12], is the 
one point in the two parameters fitting. We should consider how to select the region to 
determine the adding parameter €, How many energy points should be taken to determine 
€, In what an economic way the total luminosity should be allotted at these points to get 
the smallest statistical error of the r mass. 

Fix the total luminosity of 100pb~', in which 50pb~* is arranged to spend at the 
first point, 3.55379 GeV, to fit the 7 mass, the remainder luminosity is used for the 
€ determination. According to the distribution of the cross section, the higher energy 
region with larger cross section is more suitable to determine the efficiency. The number 
of the energy points can be studied in the region from 3.545GeV to 3.595GeV. 

To determine the number of point at the high region, m, is fixed, then fit the single 


parameter e. Figure 2 shows the variation of the Ae and S, with the adding number of 
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Figure 3: The variation of Sm, and Se with the scan of the second energy point from 3.554 
to 3.595GeV. The small boxes and the small triangles represent Sm, and Se respectively. 
The dotted line is the cross section and solid line is the derivative of the cross section. 
They have been rescaled by 0.01 and 0.001 corresponding to drawn in the same plot with 
Sm, and Se- 


points Np with two different luminosity allot precepts. The empty legends denote the 
scheme that the total luminosity is divided evenly at these points, while the solid legends 
denote the scheme that the expected event number is required to be equal at these points. 
From the one parameter m, fitting, we know that the point at the high energy region 
has few contribution to determine the m,, the variation of the statistical error Sm, just 
because of the correlation between m, and e. That is to say, the smaller the statistical 
error Se, the better the data taking scheme is. The trend of both the two schemes shows 


only one point can afford the smallest statistical fluctuation of the e. 


Now that one point is enough to determine the efficiency, a scan in high energy region 
should be done to select the point whose two statistical errors Sm, and Se are smaller. 
Figure 3 shows the distributions of the Sm, and Se when we scan this region. The cross 
section and the derivative cross section is also shown in the same plot after scaled. Clearly, 
Energy point 3.595 GeV is sufficient for the e detection. Up to now, the number of energy 
points has been confirmed as two and the optimal location should be 3.55379 GeV and 
3.595 GeV. 


The succeeding problem is how to divide the total luminosity at the two points to get 
the smallest statistical error of the r mass. Total luminosity Liot is arranged as 100 pb~! 
at the two points, namely, Lı + Ly = Lio, where Lı and Lz denote the luminosity assigned 
at the first point 3.55379 GeV and the second point 3.595 GeV, respectively. Figure 4 
shows the variation of Sm, (circle) and Se (star) with Lı. With increase the L,, the Sm, 
decreases in the beginning. That is because the first point determines the T mass, more 


luminosity reduces the error of r mass. Subsequently, the raise of Sm, at the right part of 
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figure 4 is due to the correlation of m, and e. The minimal value of Sm, is found out by 
using a curve to fit the points of figure 4. The bottom of this curve indicates that about 


75% total luminosity should be assigned at the first point. 


S, (%) 
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Figure 4: The variation of Sm, (denoted as empty circles) and S. (denoted as stars) with 


the increasing of Lı. 


Up to now, it is clear that two energy points should be taken for the two parameters 
fitting, they are Ey = 3.55379GeV, Ey = 3.595GeV. The ratio of the luminosity at the 
two points should be Lı : Ly = 3: 1. This result can be confirmed by two dimension scan 


in the energy region 3.555 ~ 3.595GeV, as shown in figure 5. 


For 75% result is on a basis of the assumption that the total luminosity is equal 
to 100pb-', we variate the value of Liot from 20pb~! to 200pb~! to check whether or 
not the luminosity assignment ratio is identical. As shown in figure 6, Every curve has 
a U shape similar with that in figure 4. The bottom of every curve reveals the best 
luminosity dividing proportion to the first point from each total luminosity Liot- The 
curve drawn below the series curves in this figure is the Sm, distribution obtained from 
the one parameter fitting. For considering the correlation between m, and e€, the same 
level of Sm, in two parameter fitting need more luminosity than that in one parameter 
fitting. Figure 7 shows the proportion of total luminosity allotted at the first point for 
each different Liot. It is a linear relationship between L; and Liot, and the slope is 0.75, 
and keeps uniform with the increase the total luminosity. So the optimal luminosity 


assignment ratio can be extended to different Liot values. 


For a summary of the two parameter, m, and c€ fitting, the optimal data taking scheme 
indicates that two points, at FE, = 3.55379GeV, Ey = 3.595GeV can acquire the smallest 
statistical error of m,. About 75% of the total luminosity should be assigned at the first 


point. 


Figure 5: The variation of S,,,.(MeV) with a 2 dimension scan of E; and E; in the energy 
region 3.55 ~ 3.595GeV. 
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Figure 6: A series curves of Sm, vs Lı for different total luminosity Liot variate form 
20pb~' to 200pb~* with a step of 10pb~'. The bottom one is the Sm, VS L curve in 1 


parameter fitting. 
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Figure 7: The optimal proportion of the luminosity allotted at the first energy point 
E = 3.55379GeV for different total luminosity Lio. 
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Figure 8: The relationship between error of tau mass and data taking point. 
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Figure 9: The variation of error of r mass and the location of energy. 


5 Three-parameter fitting 


In this section, the background cross section ogq is treated as a free parameter in 
three parameters fitting together with m, and e. The benefit of adding this parameter is 
that the correlation with other two parameters has been considered automatically. 

According to the previous T mass experiment, the background events should be a tiny 
number comparing with the observed number of events. It is acceptable to assume that 
the background cross section is a constant. We take the typical value, 0.024pb~', as the 
initial value in the simulation, which is given by reference [4]. The background region 
(from 3.50 GeV to 3.54 GeV) is supposed to contribute to the ogg determination. 

The same as two parameters fitting, we should consider how many points should be 
taken in the background region, at which energy these points should be locate, and how 
much luminosity should be arranged at these points. 

Figure 8 shows the relationship between the error of 7 mass and the number of data 
taking points in the background region. During the study the first point is at Ey, = 
3.55379GeV with Lı = 750pb~' and the second point is at Ep = 3.595GeV with Ly = 
250pb~'. The luminosity allotted on the background is 5000pb~! and is divided evenly at 
these points. According to Figure 8, increasing of the points has no effect on the accuracy 
of the 7 mass. So, one point is enough to be taken to determine the background. 

Totally, three energy points will be taken in the three parameters fitting program. 
similar as the two parameter fitting, the first point is at Ey, = 3.55379GeV with Ly = 
750pb~' and the second point is at Ep = 3.595GeV with Ly = 250pb~'. The energy of 
the third point is selected among the five points: 3.50, 3.51, 3.52, 3.53, 3.54GeV with the 
luminosity of 5000pb~'. The relationship between the error of r mass and the energy 
are shown in figure 9. It is clear that the location of the energy point in the background 
region has no impaction on the error of m, mass. So that, 3.50GeV is chosen as the third 
point. 

Then a immediate question is, how much luminosity should be arranged at the third 
point? Figure 10 shows the relationship between the error of r mass and the luminosity at 
the third point. In figure 10, the star and the dot denote the cases that the total luminosity 
at the first and the second point is 200 pb~! and 500 pb™!, respectively. There is a minimal 
value of error of r mass, after that, the increasing luminosity at the background has no 
improvement to the accuracy of r mass value. 

We fix the total luminosity as 200 pb~', and scan the luminosity allotted at the third 
point to find the optimal value. As shown in figure 11, the smallest error of r mass Sm, = 
0.0745 MeV is obtained when the luminosity equals 27.5pb~', during the scan, the ratio 
of luminosity at the first and the second point is 3 to 1. 

However, whether or not the ratio of the luminosity, 3 to 1, at the first and the second 


point is still the optimal value in the three parameter fitting? The distribution between 
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Figure 10: The relationship between error of r mass and the luminosity at the background 
point. The star line denotes that the total luminosity at the last two point is 200pb71, 
while the dotted line denotes that the total luminosity is 500pb~t. 


the error of r mass and the ratio of the luminosity at the two points is shown in Figure 
12. During the fitting, the luminosity at the third point is fixed as 27.5pb~!. As indicated 
in Figure 12, 3 to 1 is still the optimal value for the ratio of the luminosity at the first 
and the second point. 

We extend the total luminosity from 100 pb! to 500 pb~!, as shown in Figure 13, the 
optimal value of the luminosity at the third point is located in the band of about 10% of 
the total luminosity. The best ratio of the luminosity arranged at the first and the second 
point is fixed as 3 to 1. Then, fix the optimal luminosity at the third point, scan the ratio 
of the luminosity at the first and second point, the smallest errors of rT mass are obtained 
at about L,/(L; + L2) ~ 0.75, as shown in Figure 14, where L; denotes the luminosity at 
the i point. 

For a summary of 3 parameters fitting, three points are needed to determine m+, €, 
and dsg respectively. The optimal locations for the first and the second point are 3.55379 
GeV and 3.595 GeV. The best ratio of the luminosity at the two points is 3 to 1. 3.5 GeV 
is selected as the location of the third point, about 10% of the total luminosity should be 


arranged at this point. 


6 Discussion 


Figure 15 shows the comparison of the statistic error of T mass using one parameter, 
two parameters and three parameters fitting cases. Clearly, to obtain a certain statistic 
error, the three parameter fitting method need more luminosity, However the systematic 
uncertainty will be simplified. 

The designed peak luminosity of BEPCII is 1 x 10°8cm~'s~!. If the average efficiency 


of luminosity is taken as 50% of the peak value, no more than three days data taking will 
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Figure 11: The relationship between statistic error of tau mass and the luminosity at 


the background point. The total luminosity is fixed as 200 pb™t, the ratio of apportion 


luminosity at the first and the second point is 3 to 1. 
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Figure 12: The distribution of statistic error of r mass and the ratio of luminosity allotted 


at the first and the second point. 
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Figure 13: The relationship between error of tau mass and the luminosity allotted at the 
background point. The total luminosity is extended from 100 pb~! to 500 pb~!(From top 
to bottom). 
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Figure 14: The relationship between error of tau mass and the ratio of the luminosity 
allotted at the first and the second points. The total luminosity is extended from 100 
pb~! to 500 pb~!(From top to bottom). 
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Figure 15: The comparison of the statistic error of r mass using one parameter, two 


parameter and three parameter fitting methods (from left to right). 
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lead to the statistic error of r mass less than 0.1 MeV using three parameters fitting. We 
notice that this evaluation is on a basis of eu tagged events only. According to previous 
BES result, the number of multi-channel-tagged (such as ee, eu, eh, uu, wh, hh, where h 
denotes hadron) events is at least five times more that that of eu tagged events. If more 


channels are utilized to tag T-pair final states, more statics are excepted. 


7 Summary 


The data taking scheme is optimized using 2 parameters fitting of the m,, and € and 3 
parameters fitting of the m,, € and ogg. The point number and its location of the energy 
points are optimized. The precision of the statistical error of m, is less than 0.1 MeV if 
three days data taking for 3 parameters fitting. The systematic error will be studied as 


soon. 
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